A cosmological dust model with extended f(x) gravity. 
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Introducing a fundamental constant of nature with dimensions of acceleration into the theory of 
gravity makes it possible to extend gravity in a very consistent manner. At the non-relativistic 
level a MOND-like theory with a modification in the force sector is obtained, which is the limit of 
a very general metric relativistic theory of gravity. Since the mass and length scales involved in the 
dynamics of the whole universe require small accelerations of the order of Milgrom's acceleration 
constant ao, it turns out that the relativistic theory of gravity can be used to explain the expansion 
of the universe. In this work it is explained how to use that relativistic theory of gravity in such a 
way that the overall large-scale dynamics of the universe can be treated in a pure metric approach 
without the need to introduce dark matter and/or dark energy components. 
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I. INTRODUCTION 

Current cosmological data is generally explained in- 
troducing two unknown mysterious dark components, 
namely dark matter and dark energy. These ad hoc 
hypothesis represent a big cosmological paradigm, since 
they arise due to the fact that Einstein's field equations 
are forced to remain unchanged under certain observed 
astrophysical and cosmological anomalies. 

A natural alternative scenario would be to see whether 
viable cosmological solutions can be found if dark un- 
known entities are assumed non-existent. The price to 
pay with this assumption is that the field equations of 
the theory of gravity need to be extended, and so new 
Friedmann-like equations will arise. The most natural 
approach to extend gravity arises when a metric exten- 
sion f{R) is introduced into the theory [see e.g. @, and 
references therein]. 

In a series of recent articles, Bernal et al. H, HI , H ernan- 
dez, Jimenez, and Allen Hernandez et al. [12[ , Men- 
doza et al. |IJ3| have shown how relevant the introduction 
of a new fundamental physical constant ao ~ 10 _10 m/s 2 
with dimensions of acceleration is in excellent agreement 
with different phenomenology at many different astro- 
physical mass and length dimensions, from solar-system 
to extragalactic scales. The introduction of the so called 
Milgrom's acceleration constant ao in a description of 
gravity means that any gravitational field produced by 
a certain distribution of mass (and hence energy) needs 
to incorporate the acceleration ao together with Newton's 
gravitational constant G and the speed of light c in the 
description of gravity. 
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In this article we explore the consequences of the rela- 
tivistic extended metric f(x) theory of gravity [3[ applied 
to the present expansion of the universe. The article is 
organised as follows. Section Ql] gives a brief summary of 
the extended theory of gravity by Bernal et al. (H, gen- 
eralising the Newnotian description by Mendoza et al. 
[lfjj . Section IIIII interconnects this extended relativistic 
description of gravity with a metric description of grav- 
ity for which the energy-momentum tensor appears in 
the gravitational field's action. On section ITVl we use the 
developed theory of gravity for cosmological applications 
in a dust universe and see how it is a coherent repre- 
sentation of gravity at cosmological scales. Finally on 
section fVl we discuss the consequences of the developed 
approach of gravity and some of the future developments 
of the theory. 



II. RELATIVISTIC METRIC EXTENSION 

Finding a relativistic theory of gravity for which one 
of its non-relativistic limits converges to MOND yields 
usually stran ge a ssumptions and/or complicated ideas 
[see e.g. H, R Il8j. A good first approach was provided 
by a slight modification of Einstein's field equations by 
Sobouti [21|, but the attempt is not complete. 

In order to find an elegant and simple theory of gravity 
with a non-relativistic weak-field limit MONDian solu- 
tion, Bernal et al. [3[ used a correct metric interpretation 
of Hilbert's action Sf in such a way that: 



S f =- 



\QttGL\j 



/(x)v /3 5d 4 a;, 



(1) 



which slightly differs from its traditional form (see e.g. 
[ESH3) since the dimensionless Ricci scalar: 



2 



X := L 2 AI R, 



(2) 



has been introduced. In equation ([2]), R is the tradi- 
tional Ricci scalar and Lm defines a length fixed by the 
parameters of the theory. For fix) = X the standard 
Einstein-Hilbert action is obtained. The matter action 
has its usual form: 



Sty 



2c 



(3) 



with C m the matter Lagrangian density of the system. 

The null variations of the complete action, i.e. 
8 (Sh + Sm) = 0, yield the following field equations: 



fix) Xn» ~ \fix)9» v - L 2 M (V^V, - 9fW A) f(x) 



IttGL 



(4) 



M rp 



where the dimensionless Ricci tensor X/j,v is given by: 



(5) 



and R^u represents the standard Ricci tensor. The 
Laplace-Beltrami operator has been written as A :— 
V Q V Q and the prime denotes derivative with respect to 
its argument. The energy-momentum tensor is de- 
fined through the following standard relation: SS m = 
— (l/2c) T a p 5g a P . In here and in what follows, we choose 
a (+, — , — , — ) signature for the metric g^ v and use Ein- 
stein's summation convention over repeated indices. 
The trace of equation Q is: 

fix) X 2 fix) + 3L 2 M A fix) = ^T, (6) 

where T := T™. 

Bernal et al. [3( showed that the choice 



fix) = X 



- v 3/2 



(7) 



has a MOND-like solution at the weakest non-relativistic 
limit of the theory (for which x <C 1) with: 



L M = Crj.W, with C = 



2^2 



(8) 



where the "mass-length" scale Im and the gravitational 
radius r g are given by: 



1/2 



GM \ 
ao / 



GM 



(9) 



In general terms, the function fix) must satisfy the 
following limits: 



fix) 



fx 3/2 , forx<l, 
IX, 



for x> 1. 



(10) 



In other words, general relativity is recovered when x 3> 
1 in the strong field regime and the relativistic version of 
MOND with x 3 / 2 is recovered for the weak field regime 
of gravity when x <C 1. 

The mass dependence of x an d Lm means that 
Hilbert's action (JlJ is a function of the mass M. This 
is usually not assumed, since that action is thought to 
be purely a function of the geometry of space-time due 
to the presence of mass and energy sources. However, it 
was Sobouti [2l[ who first encountered this peculiarity 
in the Hilbert action when dealing with a metric gen- 
eralisation of MOND. Following the remarks by Sobouti 
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and Mendoza and Rosas-Guevara [17[ one should not 
be surprised if some of the commonly accepted notions, 
even at the fundamental level of the action, require gen- 
eralisations and re-thinking. An extended metric theory 
of gravity goes beyond the traditional general relativity 
ideas and in this way, we need to change our standard 
view of its fundamental principles. 



III. FiR,T) CONNECTION 

For the description of gravity shown in section |TT] it fol- 
lows that an adequate way of writing up the gravitational 
field's action is given by: 



Si 




16ttG 



(11) 



The function Lm is a function of the mass of the sys- 
tem and in general terms it is a function of the space- 
time coordinates. For the particular case of a spheri- 
cally symmetric space-time it coincides with the mass of 
the central object generating the gravitational field as 
expressed in equations ([8]) and ([9]). Generally speaking 
what the meaning of M would be for a particular dis- 
tribution of mass and energy needs further development, 
beyond the scope of this work. Nevertheless one expects 
that for a systems with high degree of symmetry (such 
as as a spherically symmetric spacetime or a Friedmann- 
Lemaitre- Robertson- Walker one) , the function M would 
be given by the standard mass-energy relation [see e.g. 
El: 



M := 



(12) 



The field equations produced by the null variations of 
the addition of the field's action S{ + S m can be con- 
structed in the following form. Harko et al. [Io| have 



3 



built an F(R, T) theory of gravity, so making the natu- 
ral identification: 



where the four-force 



F(R, T) := 



(13) 



J M 



it is possible to use all their results for our particular case 
expressed in equation (|13|) . For example, the null vari- 
ations of the complete action St + S m for the particular 
case of equation flT3"|) yield the following field equations 
0: 



fR ) R 1 

M / Llj M 

8ttG 



fg^v 



4 J- ,1V 



/ 



T 2 

M 



T 2 

M J T 



(14) 



A" 



8ttG 



(pc 2 + 
- u tJ -u v )V vP 
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'8ttG 






I c + 


(£>] 

\ M J T-i 



(20) 



is perpendicular to the four velocity dx"/ds. As ex- 
plained by Harko et al. [10], the motion of test particles 
is geodesic, i.e. A M = and/or V Q T Q( 3 = 0, (i) for the 
case of a pressureless p = (dust) fluid and (ii) for the 
cases in which Ft(R, T) = 0. 

In what follows we will see how all the previous ideas 
can be applied to a FLRW dust universe and so, the di- 
vergence of the energy momentum tensor in equation (|18[) 
is null. It is worth noting that this condition on the 
energy-momentum tensor for many applications needs to 
be zero, including applications to the universe at any 
epoch. 



with a trace: 



IV. COSMOLOGICAL APPLICATIONS 



f R R 2/ 



T 2 



T 2 

M 



3A 



T 2 

M 



-T — 



T 2 



r + e 

_ ' (15) 

where the subscripts R and T stand for the partial deriva- 
tives with respect to those quantities, i.e. 



BR' ^ 



For an isotropic Friedmann-Lemaitre-Robertson- 
Walker (FLRW) universe, the interval ds is given by [see 
e.g. El: 



ds 2 =c 2 dt 2 -a 2 (t) 



dr 2 



1 — nr 2 



(21) 



-2T„ 



(17) 



where a(t) is the scale factor of the universe normalised 
to unity, i.e. a(to) = 1, at the present epoch to, and 
The tensor 6^ is such that Q^Sg^ := g af) 8T a3 and for the angular displacement dfl 2 := d9 2 + sin 2 6 dip 2 for the 
the case of an ideal fluid it can be written as [I3|: po lar d9 and azimuthal d^ angular displacements with 

a comoving coordinate distance r. From now on, we as- 
sume a null space curvature k — at the present epoch in 
accordance with observations and deal with the expan- 
sion of the universe dictated by the field equations (fill) , 
avoiding any form of dark unknown component. Since we 
are interested on the compatibility of this cosmological 
model with SNIa observations, in what follows we assume 
a dust p = model for which the covariant divergence 
of the energy-momentum tensor vanishes, and so as dis- 
cussed in section [TTT] the trajectories of test particles are 
geodesic. 

For simplicity, let us consider that the function f(x) 
obeys a power-law relation: 



Note that equation (fl4|) or (fT5|) converges to the 
field (@| and trace relations as discussed in section HT1 
when one considers a point mass generating the gravi- 
tational field, i.e. when Lm = const, and so d/dR = 
L 2 M d/d X - 

In general terms, the F(R,T) theory described by 
Harko et al. [l(| produces non-geodesic motion of test 
particles since: 



~ \ I 2 



J M / T 



T 2 

M / T 



fix) = x h 



(22) 



/ 



(T^ + e^JVIn ) (i?,T) + V^6^ 



J M / T 



(18) 



and as such the geodesic equation has a force term: 



with a constant exponent b. Let us now rewrite the field 
equations (I14[) inspired by the approach first introduced 
by Capozziello and Fang Q (see also Capozziello and 
Faraoni M]) as follows: 



d 2 x>" 



ds 2 



(19) 



8ttG 



8vrG T F R 



rriCUrV 



(23) 
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where the Einstein tensor is given by its usual form: 



(24) 



and 



T-icurv 



— { -(F-RF R )-AF R ) 9 ^ 



(25) 



represents the "energy-momentum" curvature tensor. 
Since Too = pc 2 , then it will be useful the identifica- 
tion Too := PcurvC 2 . With this last definition and using 
the fact that the Laplace-Beltrami operator applied to a 
scalar field -0 is given by [see e.g. [3]: 



1 



V^g 



(26) 



LOC p c 



and 



(F - RF R ) /2 + d 2 F R /dt 2 + 3HdF R /dt 
c 2 (RF R — F) /2 — 3HdF R /dt - 



(32) 



(33) 



On the other hand, note that the mass M that appears 
on the length Lm must be the causally connected mass at 
a certain cosmic time t, since particles beyond Hubble's 
(or particle) horizon with respect to a given fundamental 
observer do not have any gravitational influence on him. 
In other words: 



,1/ = 



4 (f 

-9 d x = ^Pjp' (3 4 ) 



where 



then 



8irGF R 



1 (up p\ * HdF X 
~(RF R -F)-^ — 



(27) 



where H := a(t)/a(t) represents Hubble's constant. 

With the above definitions and using the 00 component 
of the field's equations (j2"3")l and the relation [cf. 



R 



KC 
~^2~ 



H(t) 



(35) 



is the Hubble radius or the distance of causal contact at 
a particular cosmic epoch [l5[ . From this last relation it 
follows that the length ([5]) is given by: 



ca 



1/4 H 9/i ■ 



(28) and so, by using relation (|22l) then: 



(36) 



between Ricci's scalar and the derivatives of the scale fac- 
tor for a FLRW universe, then the dynamical Friedman's- 
like equation for a dust flat universe is: 



IttG 



— — Ft) — 

8ttG T ) F r 



Pcurv 



(29) 



The energy conservation equation is given by the null 
divergence of the energy-momentum tensor: 



f 8nG 
\ c 4 



F T )(p+ SHp) 



-P- 



dF T 
dT' 



(30) 



For completeness, we write down the correspondent 
generalisation of Raychadhuri's equation for a dust flat 
universe: 



a 
2- 

a 



inGp c 



(31) 



where the "curvature-pressure" 



dF K 
dt 



= b{b-l)R"- L L 



6_l r 2(6-i; 



M 



H 



j-q-2 3 
1 -q 2 



dFT - 3 (h 

dt 2 { ' pc 2 ' 



where 



a 
(37) 

(38) 



, . 1 d 2 a o , 1 d 3 a o , . 

*® l= -a& H ' &nd ,:= ^ r ' (39) 

are the deceleration parameter and the jerk respectively. 

As it is usually done, let us find power law solutions 
satisfying 



a(t) = o(t ) ( - 



P(t) = Po 



a(t ) 



(40) 



for the unknown indices a and /3. With these and the 
value of Lm from equation (|36|) . the curvature den- 
sity (|57]) is 



5 



3H 2 
8nG 



(6-1) 



(1 - 9) -^^-2( v /3+ « 



( 41 ) 

Substitution of the previous relations on Friedmann's 
equation (p9|) it follows that: 



H 



, _ 87TG/9 



(42) 



where 



Z:=l + (&-!) 



J-g-2 4(l-g) ,3/3 
1 - g 6 2 \ ex 



(43) 

is a dimensionless function. 

An important result can be obtained evaluating equa- 
tion (|4"2"j) at the present epoch, yielding: 



an 



(36_ 5 )/(b_l) 



"(44) 

where the density parameter O^tt a ^ the present epoch 
has been defined by it's usual relation: 



n 



(0) 



3H 2 p 
8ttG ' 



(45) 



In other words, the value of Milgrom's acceleration 
constant a at the current cosmic epoch is such that 



a « c x ff . 



(46) 



The numerical coincidence between the value of Mil- 
grom's acceleration constant ao and the multiplication 
of the speed of light c by the current value of Hubble 's 
constant Hq has been noted since the early development 
of MOND [see e.g. [9j, and references therein]. Note that 
equation (|46p means that this coincidence relation occurs 
at approximately the present cosmic epoch in complete 
agreement with the results by Bernal et al. [2J where it is 
shown that ao shows no cosmological evolution and hence 
it can be postulated as a fundamental constant of nature. 

For the power law l|22p and the assumptions made 
above, it follows that the energy conservation equa- 
tion (|30|) is given by: 



(p + 3Fp) 
where: 



— ( A- 

BttG V 9 



B H) R b Lf- 1] = 0, (47) 



B :=- 



9 6_ 
2 



1 27(6-1) 



3 6(6 
2 



i)(i 



Direct substitution of the density power law (|40|) into 
relation (|47[) gives a constraint equation between a, /3 
and 6: 



0=- 



56 



a V 36 — 5 



(48) 



Let us now proceed to fix the so far unknown param- 
eters of the theory a, f3 and 6. To do so, we need reli- 
able observational data and as such, we use the redshift- 
magnitude SNIa data obtained by Riess et al. [20( and the 
following well known standard cosmological relations [see 
e.g. EH: 

1 + z = a(t )/a(t), (49) 
H (z) = 5 log 10 [H d L (z)] - 5 log 10 h + 42.38, (50) 



d L (z) = (l + z) 



o H(z) 



dz, 



(51) 



for the cosmological redshift z, the distance modulus fi, 
the Luminosity distance and where the normalised 
Hubble constant h at the present epoch is given by 
h := H Q / (l00kms~V M P c ). A l so > lrom equation (gDJ it 
follows that 



H{a) = H Q 



a(t ) 



-l/c 



Hoil + z) 1 ^, 



(52) 



and the substitution of this into equation (1511) gives the 
distance modulus d^ as a function of the redshift z. This 
means that the redshift magnitude relation (|50j) is a func- 
tion that depends on the values of the current Hubble 
constant Hq and the value of a. Figure [1] shows the 
best fit to the redshift magnitude relation of SNIa ob- 
served by Riess et al. [20(, yielding a = 1.359 ±0.139 and 
h = 0.64± 0.009. The best fit presented on the figure was 
obtained using the Marquardt-Levenberg fit provided by 
gnuplot (http://www.gnuplot.info) for non-linear func- 
tions. This value does not conclude the whole description 
of the problem, since /3 and 6 are still unknown but ac- 
cording to the constraint equation (|48p only one of them 
is needed in order to know the other once a is known. 

The parameter /3 can be found from conservation of 
mass arguments, since the total mass of the universe 
M to t = 47r Jq " 1 ^ pr 2 a(t) 3 dr = const., where the upper 
limit of the integral is the radius of the whole universe. 
Since a(t) and p(t) are time dependent functions, the only 
way the mass of the universe is conserved is by requiring 
a 3 p = const, and so, j3 = 3. This argument is exactly the 
one used in standard cosmology when dealing with a dust 
FLRW universe [see e.g. [l5|. Using this value of j3 and 
the one already found for a, it follows that 6 = 1.57±0.56 
which is within the expected value of 6 = 3/2 discussed 
in section |TTJ 

For completeness, we write down a few of the cosmog- 
raphycal parameters obtained by this fix) gravity ap- 
plied to the universe: 
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FIG. 1. Redshift magnitude plot for SNIa showing the dis- 
tance modulus fi as a function of the redshift z for SNIa as 
presented by Riess et al. [23|. The dotted red line shows the 
best fit to the data with the f(x) gravity theory applied to a 
flat dust FLRW universe (see text) with no dark components. 
The continuous blue line represents the best fit according to 
the standard concordance dust ACDM model. 



h = 0.64 ±0.009, q = -0.2642 ±0.075, 
jo = -0.1246 ±0.004. 



(53) 



characteristic size r is such that x := Im/t ^5 1 i s i n 
the MONDian regime of gravity. For the case of the 
universe, x ~ a few, and as such if not totally in the 
MONDian regime of gravity, then it is far away from 
the regime of Newtonian gravity. The relativistic version 
of this means that the universe is close to the regime 
for which fix) — X 3 ^ 2 an( l so b — 3/2. This is a very 
important result since, seen in this way, the accelerated 
expansion of the universe is due to an extended gravity 
theory deviating from general relativity. It is quite inter- 
esting to note that the function fix) = X 3 ^ 2 which at its 
non-relativistic limit is capable of predicting the correct 
dynamical behaviour of many astrophysical phenomena, 
is also able to explain the behaviour of the current accel- 
erated expansion of the universe. 

Seen in this way, the behaviour of gravity towards the 
past (for sufficiently large redshifts z) will differ from 
fix) = X 3 ^ 2 ancl eventually converge to f(x) = X: i- e - 
the gravitational regime of gravity is general relativity for 
sufficiently large redshifts. A very detailed investigation 
into this needs to be done at different levels in order to be 
coherent with many different cosmological observations 
[see e.g. [HI . 

It is quite remarkable that a metric extended theory of 
gravity is able to reproduce phenomena from mass and 
length scales associated to the solar system up to cosmo- 
logical scales. Many other cosmological applications of 
the theory will be addressed elsewhere. 



DISCUSSION 
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